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We propose a physical mechanism to explain the crystal symmetries found in macromolecular and 
supramolecular micellar materials. We argue that the packing entropy of the hard micellar cores is 
frustrated by the entropic interaction of their brush-like coronas. The latter interaction is treated as 
a surface effect between neighboring Voronoi cells. The observed crystal structures correspond to the 
Kelvin and Weaire-Phelan minimal foams. We show that these structures are stable for reasonable 
areal entropy densities. 



PACS numbers: 83.70.Hq 82.70.-y 61.50.Ah 

Dendritic polymers fl]-^, hyper-branched star poly- 
mers |^,|| and diblock copolymers (6|j^] represent a new 
class of molecular assemblies all of which form a vari- 
ety of crystalline lattices, many of which are not close- 
packed. These assemblies are all characterized by com- 
pact cores and brush-like, soft coronas. These systems 
might be modeled by treating the micelles as sterically 
interacting hard spheres and it would follow that their 
crystalline phases should be stackings of hexagonal-close- 
packed (HCP) layers. Recently J|j it has been shown that 
the face-centered cubic (FCC) lattice maximizes the to- 
tal entropy and so hard-sphere crystals should form FCC 
structures. Note that the entropic difference between the 
various HCP lattices is a global issue: the local arrange- 
ment of spheres is the same for all close-packed variants 
and thus the lattice cannot be predicted from nearest- 
neighbor interactions. In order to understand the rich- 
ness of crystal symmetries in the micellar systems, we 
propose an additional global consideration: we add an 
interaction proportional to the intcrfacial surface area 
between the cages which contain each micelle (Voronoi 
cells). Though approaches based on self-consistent field 
theory and two-body interactions can yield non-close- 
packed lattices , we propose a universal explanation 
for a host of new structures and present a new paradigm 
for the rational design and control of macromolecular as- 
semblies |IT| . 

The interfacial interaction arises through the entropy 
of the brush-like coronas of the micelles. Because of 
constraints on their conformations, the brushes suffer an 
entropic penalty proportional to the interfacial area be- 
tween the Voronoi cells surrounding each sphere. Thus 
they favor area-minimizing structures, precisely the type 
of structures that dry foams might make. Over a cen- 
tury ago, Lord Kelvin proposed that a body-centered- 
cubic (BCC) foam structure had the smallest surface-to- 
volume ratio |l2| but in 1994 Weaire and Phelan found 
that a structure based on the A15 lattice jl3| was more 
efficient. We note that neither the BCC nor A15 struc- 
tures are close-packed and thus there is a fundamental 
frustration between the hard-core volume interaction and 
the surface interaction due to overlapping soft coronas. 

For concreteness, in this paper we focus on structures 
observed in a family of dendrimer compounds consist- 



ing of a compact poly(benzyl ether) core segment and 
a diffuse dodecyl corona jj]]^]. These conical dendrimers 
self- assemble in spherical micelles which subsequently ar- 
range into the A15 lattice (Fig. Q). The interaction be- 
tween the micelles is primarily steric, i.e., repulsive and 
short-range. The micellar architecture suggests that the 
potential is characterized by three regimes. At large dis- 
tances, the micelles do not overlap and the interaction 
vanishes. As the coronas begin to overlap, the entropy 
of the brush-like coronas decreases, which gives rise to 
a soft repulsion between the micelles. Finally, at small 
separations the coronas begin to penetrate the compact 
cores: this is very unfavorable and gives rise to hard-core 
repulsion. This energy landscape is in qualitative agree- 
ment with recent, detailed molecular dynamics simula- 
tions in. 
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FIG. 1. Various lattices: a) Face-centered cubic, b) 
body-centered cubic, c) A15 lattice, and d) columnar rep- 
resentation of A15 lattice. In the A15 lattice, columnar and 
interstitial sites are drawn in grey and black, respectively. 

Although both originate in steric interaction, the two 
repulsive regimes are characterized by very different func- 
tional behaviors. The hard part of the potential results 
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in a restricted positional entropy of the micelles which 
depends on the free volume, the difference between the 
actual and the hard-core volumes. The soft part comes 
from the decreased orientational entropy of the chains 
within the overlapping coronas. The matrix of overlap- 
ping coronas can be thought of as a compressed bilayer 
and thus the free volume may be written as a product of 
the interfacial area A and the average spacing between 
the hard cores d so that at any given density 



Ad — constant. 



(1) 



Though this approximation ignores the curvature of 
brush-like coronas, the dendrimers are relatively close 
and we expect this constraint to hold in this system. 
Since the repulsion decreases monotonically with dis- 
tance, the system will favor a maximum thickness d and 
will thus tend to minimize the interfacial area. Hence 
our proposed interfacial interaction, which is incompati- 
ble with the bulk free energy minimized by a close-packed 
arrangement of micelles. In the following, we compare 
the free energies of FCC, BCC, and A15 lattices and 
estimate the strength of the interfacial interaction such 
that the structure of the micellar crystal is dictated by 
the minimal-area principle. 

The calculation of the bulk free energies of condensed 
systems is fairly complicated even for hard-sphere sys- 
tems and the best theoretical results are obtained nu- 
merically. It is interesting to note that elaborate ana- 
lytic models, such as the high-density analog of the virial 
expansion jlql and the weighted-density-functional ap- 
proximation fllql , are only slightly better than the simple 
cellular free- volume theory |l6],[l7| . The free- volume the- 
ory is a high-density approximation where each micelle is 
contained in a cell formed by its neighbors, and the com- 
munal entropy associated with the correlated motion of 
micelles is neglected. 

Within this theory, the positional entropy of a micelle 
is determined by the configurational space of its Voronoi 
or Wigner-Seitz cell. In the FCC lattice, the centers of 
mass of the micelles are within rhombic dodecahedra [18l , 
while in the BCC lattice they are contained in regular 
octahedra although the BCC Voronoi cell is an orthic 
tetrakaidecahedron For these lattices, the bulk free 
energy of a micelle is given by 



bulk 
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where X is either FCC or BCC, n — pR 3 is the reduced 
number density, and R is the hard-core radius of micelles. 
The coefficients a FCC = 2 5 / 2 and a BCC = 2 2 3 1 / 2 depend 
on the shape of the cells, whereas /3 FCC = 2~ 5 / 6 and 
/3 BCC = 2 -5 / 3 3 1 / 2 are determined by their size. 

The A15 lattice is somewhat more complicated: as 
shown in Fig. |l|d the A15 unit cell includes 6 colum- 
nar sites, which make up 3 perpendicular interlocking 
columns, and 2 interstitial sites. A pseudo- Voronoi 



construction (subject to the constraint that all cells 
have equal volume) for this lattice leads to a parti- 
tion consisting of irregular pentagonal dodecahedra and 
tetrakaidecahedra with two hexagonal and twelve pentag- 
onal faces |nj . Because of the irregularity of the cells, we 
calculate the bulk entropic free energy numerically, and 
the result is shown in Fig. ^. For our purposes we re- 
quire an analytic form: by substituting the dodecahedra 
and tetrakaiecahedra by spheres and cylinders, respec- 
tively, and allowing for two adjustable parameters C and 
S, which measure the deviation of the Voronoi cells from 
these spheres and cylinders, we find 



Ft 



A15_ 



bulk- 



-k R T 



1 (AttS ( a/5 
4 \ 3 v 4nV3 



2nV3 



) 



(3) 



This form is within 0.1% of the numerical result with 
S w 1.64 and C w 1.38. 
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FIG. 2. Bulk free energies of FCC, BCC, and A15 lattice 
at reduced densities above the melting point. Solid lines: an- 
alytical results, Eqs. (^) and (^|); circles: numerical results. 

The interfacial free energy is minimized by the division 
of space with smallest area. The problem of finding the 
partition of space into equal-volume cells with the mini- 
mum interfacial area was first studied by Kelvin ]2C| ] : he 
proposed a BCC lattice of orthic tetrakaidecahedra with 
slightly curved hexagonal faces to satisfy the Plateau 
rules 0|. However, the Weaire-Phelan partition, which 
differs from the equal-volume Voronoi construction for 
the A15 lattice only in a delicate curvature of the pentag- 
onal faces, is 0.3% more efficient We note that the 
BCC and A15 structures are among the simplest tetrago- 
nal close-packed lattices [^l], suggesting that other, more 
complex close-packed lattices might be more efficient 
still. However, the A15 structure appears to be the most 
efficient, although no proof of its supremacy exists. 

To argue that the A15 and BCC lattices are the equilib- 
rium structures formed by micelles, we must estimate the 
entropy penalty per unit area and translate this into an 
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entropy per dodecyl chain. The dodecyl bilayer is mod- 
eled as a polymer brush consisting of chain molecules at- 
tached to hard cores, and in the limit of high interdigita- 
tion its free energy consists solely of the excluded- volume 
repulsion of the chains: 



-^surf — 
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where d is the layer thickness, I is a parameter with the 
dimension of length and N is the number of chains per 
micelle p3]. Since the bilayer must fill the free volume, 
Am<1 = 2(n~ 1 — 47r/3)-R 3 , where Am is the interfacial 
area per micelle. Thus the interfacial free energy of a 
micelle is 
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where -/ cc = 2 5 / 6 3 = 5.345, 7 BCC = 5.306, and 7 A15 = 
5.288 [g§. 

We now calculate the range of £ such that the total 
free energy 
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is minimized by the BCC and A15 lattices rather than 
by the naive, close-packed, FCC lattice. In order to es- 
timate the strength of the soft repulsion, we must first 
determine the actual reduced density n. Since the hard- 
core radius of the micelles is unknown, we limit n by 
recognizing that it must be larger than the melting den- 
sity, n « 0.120 for hard spheres |lq| , and that it must be 
smaller than the close-packing density of the A15 lattice, 
n = 0.125. The most conservative lower bound of £ cor- 
responds to the lowest possible density, i.e., the melting 
density. With N = 162 chains per micelle Q, we find 
that at n = 0.120 the FCC to BCC transition occurs 
for £ w OAR and the BCC to A15 transition occurs for 
£ w 0.3i?. This corresponds to an entropy per chain of 
0.5k B and 1.5k B , respectively. Both values are of the cor- 
rect order of magnitude and the higher value of the latter 
is consistent with the relative rarity of the A15 phase. 

Since we expect that each chain has at least k B of en- 
tropy, we conclude that the energetics of the dendrimer 
micelles is dominated by interfacial effects. This is hardly 
surprising. The number of degrees of freedom of each mi- 
celle is quite large and the bulk free energy only depends 
on the position of the micelle as a whole. Since the mi- 
celles are soft, the internal degrees of freedom such as the 
chain conformations play an important role. 

This paradigm - which shows that the minimal surface 
problem can be fruitfully transplanted to the microscopic 
level - explains the morphology of a number of dense mi- 
cellar systems. The same ideas can be applied to poly- 
meric micelles made of, e.g., polystyrene-polyisoprene di- 
block copolymers dispersed in decane df]]. ^ n this case, 
the micelles are characterized by highly concentrated 
polystyrene core and diffuse polyisoprene corona, and 



they form BCC or FCC lattices, depending on the rel- 
ative length of the polystyrene and polyisoprene chains. 
The BCC lattice is observed in diblock copolymers with 
similar lengths of core and coronal segments, whereas the 
FCC lattice occurs whenever the corona is thin compared 
to core. This is consistent with our model of the impen- 
etrable core which is responsible for the hard part of the 
repulsion and which favors arrangements with large free 
volume. In addition, our model suggests that the A15 
lattice is the ground state of an asymmetric diblock with 
an exceptionally large corona or, equivalently, a corona 
made of very floppy, "entropy-rich" chains. We note that 
distinguishing between A15 and BCC in powder-averaged 
diffraction is delicate: the first three BCC reflections |24] 
are at v2> v4, and v6, while the first four A15 reflec- 
tions are at v2j v4, \/5, and and thus a careful 
study would be necessary. 

The existence of the A15 lattice in the dendrimer ag- 
gregate also may be regarded as an experimental verifi- 
cation of the recent theoretical developments in minimal 
surfaces and, in particular, Weairc and Phelan's conjec- 
ture that this structure solves the Kelvin problem. At 
this juncture, the presumably ideal A15 structure has not 
been observed unambiguously on a macroscopic scale in a 
soap froth |25}| . Last but not least, let us note that similar 
structures have been found in lyotropic materials, e.g., in 
lipid bilayers in water j2f|. In such systems the intermi- 
cellar potential also results in an effective interfacial free 
energy although it is not steric but substance-specific, 
and thus transcends the scope of this discussion. 

Our model may be further refined by including the 
effects of the curvature of the brush-like coats, the strain 
of the coronas into the interstitial regions, and solvent 
effects. In addition, the dual problem of determining the 
structure of foams might be amenable to our analysis 
through the introduction of excluded volume interactions 
between the bubbles |27| ]. Recent work focusing on two- 
body interactions has shown that 

We hope that this study elucidates the relation be- 
tween interaction and structure in supramolecular assem- 
blies. By including an additional global contribution to 
the free energy we provide a rough yet universal quanti- 
tative guideline for the design of self-organized soft mate- 
rials, which can be used for a number of applications such 
as photonic band-gap materials Q , Bragg switches [^9| , 
and porous microreactors p0[ |. By tuning the ranges of 
hard and soft repulsion, one should be able to choose 
among the spectrum of symmetries from the lattice with 
minimal interfacial area to the lattice with maximal pack- 
ing fraction and engineer the crystal structure most fitted 
for a particular application. 
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